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Abstract 



o ■ We consider a Schwarzschild black hole immersed in a ghost condensate background. It is shown 

^ that the Hawking radiation in the quanta of small perturbations around this background is highly 

suppressed- in particular it is not given by a thermal spectrum. This result is in accord with obser- 
vations that such black holes can be used to violate the generalized second law of thermodynamics, 
(^ . and thus cannot have a standard entropy /area relation. 

(N 

> 

m 

in 

p 

O 
00 
O 

> 

X 



* Electronic address: lbfeldste@bu.edul 



I. INTRODUCTION 

The classic results of Bekenstein [l| and Hawking [2| that black holes have an entropy 
proportional to their area and that they radiate as blackbodies, have proven to be of fun- 
damental importance to our current understanding of quantum gravity. It is now widely 
relieved that any complete theory of quantum gravity must satisfy the holographic principle 
3|, IJ] , which states that the number of fundamental degrees of freedom contained within a 
given region scales as the area bounding that region, not as its volume. The holographic 
conjecture itself was motivated by the noted discoveries in black hole thermodynamics, and 
has gained further support from the AdS/CFT correspondence [5|]. 

It is very interesting in this context to try to find theories for which black hole thermo- 
dynamics does not behave in the way expected from holography. If we can find low energy 
effective theories with this property, then we may be able to reach the striking conclusion 

that these theories are incompatible with fundamental principles of quantum gravity. 

i' 

An important example can be found in [6], in which the authors studied a Schwarzschild 
black hole immersed in a ghost condensate background [3]. In the ghost condensate, a scalar 
field has a non-vanishing expectation value for its kinetic term, while having essentially 
zero energy-momentum. This background breaks Lorentz invariance, and adding couplings 
between the ghost condensate and additional fields, we can obtain particles with maximum 
speeds which differ from the speed of light. 

In the black hole background, such particles have horizons at different radii than the 
usual horizon, and by extension. Hawking spectra with temperatures different than the 
usual temperature. Of course, already the thermodynamic picture of such a black hole is 
a bit mysterious- it is not clear for example what temperature (or entropy) is the "correct" 
one. Indeed, the authors of [6] showed that it is possible to use such a black hole and 
an arrangement of shells at various temperatures to violate the generalized second law of 
thermodynamics ( 8| used an alternative mechanism to reach a similar conclusion)^. 



^ This author is somewhat concerned that the perpetual motion machine discussed in [6(] (and also the one 
in [8[) has a possible loophole, in that the flux of particles heading into the black hole may cause a relevant 
shift in the ghost condensate background. This could lead to an extra flux of ghost condensate energy 
into the black hole, and therefore also to growth in the size of the horizon. This effect is at the right order 
of magnitude to possibly (but not necessarily) ruin the claimed GSL violation. To analyze this issue at a 
greater level of precision will require further study. 



These facts could lead one to speculate that the ghost condensate could not emerge in 
a low energy effective field theory coming from a consistent quantum theory of gravity. In 
fact, a UV completion of any kind for the ghost condensate has yet to be found. 

It is then important to understand what exactly the fundamental property of the ghost 
condensate theory is which causes its apparent difficulties with quantum gravity. The likely 
candidate for such a property is the violation of the null energy condition present in these 
theories.^ It is this feature which allows the ghost condensate to break Lorentz invariance 
while simultaneously carrying negligible energy-momentum- the key feature which allows 
the construction of the perpetual motion machines of [6| and Sj. 

It is interesting that, independent of any considerations involving the null energy condi- 
tion, the arguments of |6| seem to demonstrate a pathology with black hole thermodynamics 
in the ghost condensate using only the properties of the Hawking radiation in various species. 
It seems extremely interesting that one could compute a Hawking spectrum and conclude 
that a particular low energy effective field theory is incompatible with principles of quantum 
gravity. 

With this motivation, in this paper we will attempt to complete the story of the Hawking 
radiation in the ghost condensate by computing the spectrum of the ghost condensate quanta 
themselves. Indeed, this will allow us to see a pathology in the Hawking spectrum without 
having to require specific sorts of couplings to additional fields: We will find that the 
Hawking radiation in these particles is highly suppressed, and in particular, is non-thermal. 
The radiation vanishes in the limit that 217 ^ ^ ^ ^pi' where M is the black hole mass, 
and A is the UV cutoff of the ghost condensate effective theory. Moreover, the flux in 
a frequency u vanishes in this limit even if we hold flxed the product uM, thus ruling 
out a typical Boltzmann distribution such as e~'^^^, with T oc 1/M. Note that here and 
throughout the paper, we will work in natural units with Mpi = 1. 

In section ([TT]), we will review a useful calculation of the Hawking spectrum for a massless 
scalar in 2 dimensions, while in section (lllip we will review relevant aspects of ghost con- 
densation. In section ( JIVI) we will discuss the possible choices of higher derivative operators 
in the theory, and their possible impact on our results. In section (JV]) we will calculate the 



^ The null energy condition is the requirement that, for any future pointing null vector e^, T^^e^e'^ ^ 0, 
where T^v is the energy-momentum tensor. Physically this can be thought of as the requirement that 
light rays are always caused to be focused, rather than anti-focused, by sources of energy. 



Hawking spectrum of ghost condensate quanta, showing that it is non-thermaL We will 
conclude in section fIVll). 



II. MASSLESS SCALAR HAWKING SPECTRUM IN 2D 

Here we review a particular method for calculating the Hawking spectrum for a massless 
scalar field in 1+1 dimensions |9]. This method will later be shown to be generalizable to 
a calculation of the spectrum of ghost condensate quanta. Much of our analysis follows the 
line of argument in the excellent review by Jacobson [10]. 

We consider a massless scalar field ip propagating in a two-dimensional Schwarzschild 
black hole background. The metric is 

ds'=[l-—jdt'-^^—^dr'. (1) 

We make a change of time coordinate to 




SO that the metric takes the form 

ds^ = (l - v^) dr^ + 2v drd r - dr^ (3) 

where v{r) = —\/^- The velocity vector of an observer freely falling into the black hole 
from rest at infinity is always perpendicular to the surfaces of constant r. In fact r measures 
the proper time for such observers, and the (r, r) coordinate system is perfectly well behaved 
as one crosses the (future) event horizon. This coordinate system is thus well suited to a 
calculation of the Hawking spectrum generated at the horizon, and in fact it is particularly 
convenient for two additional reasons: First, the time r does not appear explicitly in the 
metric, and second, as will be important later, the ghost condensate background in the 
presence of a black hole satisfies </? oc r. 
The action for the massless scalar is 

S = j ^/^g^'V.i^V^i^. (4) 

In the metric (3) this takes the form 

S= fdrdr [{{dr + vdrW - (dr^/j)^] , (5) 



and the equation of motion 

g^^'V^V.ij = (6) 



takes the form 

dl + 2vdrdr - (1 - V^)dl + —dr + 2t; — 9, 



^ = 0. (7) 



dr dr 

As a result of the metric being independent of r, there exist solutions to the equation 
of motion of the form e~*'^'^/(r)'^. At large distances from the hole, we may take /(r) 
to approach~ e^^*^"", corresponding to outgoing or ingoing modes. A particle seen by an 
observer at infinity may be taken to be a wave packet formed out of such solutions peaked 
around some particular frequency of interest. Note that purely outgoing modes will also have 
a boundary condition that they vanish behind the horizon; without this boundary condition, 
the solution would also describe a separate bit of matter falling into the singularity. We may 
then ask, for example, about the expectation value of the number operator corresponding 
to outgoing modes of this type, with a non-zero result indicating particle production. 

Now, associated with the wave equation (Q, we have an inner product defined on complex 
solutions given by 

< si, S2>=i I n^'[slV^S2 - saV^Si]. (8) 

Here S is a given space-like hyper-surface, with n^ being its unit normal vector. Note that 
this inner product is not positive definite, although it is independent of the surface chosen 
due to the Klein-Gordon equation (l6l). On a surface of constant r, the inner product takes 
the form 

< Si, S2 >= i / dT[sl{dr + Vdr)s2 - S2{dr + f(9,.)s*]. (9) 

Now, given a solution to the equation of motion s we can define an associated operator by 

a(s)=<s,^>, (10) 

with a Hermitian conjugate 

at(s) = - <s*,^> . (11) 

Here \1/ is the quantum operator corresponding to ^. It then follows from the canonical 
commutation relations that 

[a{si),a\s2)]=<Si,S2> . (12) 



By equation (2), these are also solutions of fixed killing frequency. 



Thus if s is a solution with positive norm in the inner product ([8]) (and we normahze it 
appropriately), then the operator a{s) is in fact an annihilation operator associated with 
the given mode. If s has negative norm, then its complex conjugate has positive norm and 
will be associated with an annihilation operator. 

To compute the expectation value of the number operator associated with a given solution 
s, we must thus evaluate < a^{s)a{s) >. Note however that since the inner product 
(IH]) is conserved, we are free to perform the computation on any space-like slice of our 
choosing. Exploiting this freedom, we may thus propagate an outgoing wave packet of 
interest backwards in time until it is highly blue shifted and squeezed up against the horizon. 
The reason this is fruitful is that, while we would claim to have no a priori knowledge of 
the particle content of low energy modes (such as the one we started with), we can say 
something about the content of high energy modes: Since the evolution of the black hole 
background is taking place on timescales of order M, modes with energy much greater than 
1/M are expected to be in their ground states [ll|]. Note that this argument applies to the 
energy as seen by observers who free fall into the black hole from infinity. Equivalently, 
it applies so long as the spacial slices we are using to define the time evolution and energy 
have an extrinsic curvature which is always < 1/M. 

We may thus propagate a low energy wave-packet of interest back in time until it is 
squeezed against the horizon and composed primarily of high energy modes as seen by a 
free falling observer. Note that such high energy modes effectively propagate in flat space 
since their wavelengths are much less than the Schwarzschild radius. In fiat space, we have 
plane wave modes of the form e~*'^*e^*'^^, and in the inner product ([H]) these have positive 
norm for u > and negative norm for u < 0. When we evaluate < a^(s)a(s) >on an early 
time slice, we will thus decompose s into positive and negative free fall frequency pieces s^ 
and s_, with positive and negative norm. We will then have 

a(s) =< s, ^ >=< s_^, ^ > + < s_, ^ >= a{s^) - a"^(sl), (13) 

where in the last equality the positive and negative norms of s+ and s_ imply that a(s^_)and 
a^(sl) are bona fide annihilation and creation operators, respectively.'* 



"* More generally, there is also a component at early times which is a positive norm, low energy, ingoing 
mode, having been reflected off the geometry back out to infinity. For the wave-packets considered in 
this paper the reflected pieces will be negligible. 



This allows us to calculate the expectation value of the number operator we are interested 
in; 

< a\s)a{s) > = < a{s*_)a\s*_) > 
= < s*__,s*_ > 
= - <s„,s_ >, (14) 

where in the first line we have used the fact that the annihilation operators for the high 
energy modes annihilate the vacuum, and in the second we have used the commutation 
relation flT2|) . Thus the number operator for a low energy outgoing wave-packet is given by 
the norm of its negative free fall frequency component at early times. To find this negative 
free fall frequency part, one may explicitly take components on a basis of modes using the 
inner product (JHl). It is also possible to use a rather nice trick due to Unruh [12] as we will 
shortly review. 

As mentioned above, we would like to form our wave-packet out of solutions to the 
equation of motion of the form e~*'^'^/(r), and which become outgoing plane waves at large 
distances. Now, it is important to realize that although the coordinate r does indeed 
measure the proper time for free falling observers, modes of this form do not have definite 
free fall frequency. This is a result of the existence of the off diagonal term in the metric 
([3]); derivatives in the direction of a free falling observer's world-line are actually given by 
dr + vdr (by contrast, dr does take a derivative in a direction perpendicular to this). Of 
course, a function of the form e~"^'^/(r) is not generally an eigenfunction of this derivative 
operator. 

In flat space, our massless scalar satisfies the dispersion relation u"^ = k"^. In the curved 
space-time of the black hole, we might hope to find solutions to the equation of motion by 
using a WKB approximation. Such solutions take the form 

s^(r) oce-^^V^' '=("') ^'■', (15) 

with k{r) being given by the local solution to the dispersion relation at the given radius. 
That is, since {dr + vdr)suj = —i{uj — vk)si^, the "local" free fall frequency is given by 

uJs{r) = LU — v{r)k{r), (16) 

and we take 

Mrf = k{rf. (17) 



These WKB solutions will be good approximate solutions to the equation of motion (jTj) so 
long as the distance scales over which the geometry and the wavelength change are much 
longer than the wavelength itself. This corresponds to the requirements 

v'{r) < \k{r)v{r)\, (18) 

and 

|A;'(r)| <A;(r)2. (19) 

In the case of the massless scalar, fITBl) and flTTl) yield 

Hr) = — ^, (20) 

1 + f (r) 

where we have chosen the solution which becomes an outgoing plane wave at infinity. In 
fact for the massless scalar field, equations flTSl) and flT9|) are not satisfied unless uj ^ ^/M, 
and one would expect the WKB solutions to be invalid at wavelengths of interest of order 
the Schwarzschild radius and larger. It turns out, however, that solutions of the form flTSl) . 
with k{r) given by (J20l) just so happen to be exact solutions in the case of the massless scalar 
in 1+1 dimensions (the terms in the equation of motion one would want to drop by the 
WKB approximation actually happen to cancel each other). The WKB method will still 
be useful in the case of the ghost condensate, since in that case the analogues of equations 
(|T8|) and ^ will be valid. 

We now form a wave-packet s out of the WKB modes, and propagate it back in time 
until it is squeezed up against the horizon. We then wish to find the positive and negative 
free fall frequency pieces. To this end we thus only need to know the behavior of our modes 
close to the horizon, and we expand r = 2M + x, finding at leading order 

k ^ '-M^, (21) 

X 

and 

The constant overall amplitude is irrelevant and has been dropped. These expressions are 
valid for x > 0; recall that behind the horizon the outgoing modes we are interested in vanish 
by definition. Now, an observer falling freely into the black hole from infinity follows a path 
with r = (|A/2M(i? — t))^/^ for some constant R. Close to the horizon these paths have 
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X ~ — Ar, and thus the leading behavior for the mode (122|) along the free falling observer's 
world- line is given by s^^ ~ e*4Af(^Log[-AT] ^^^ ^^ negative, with s^j = for Ar positive. 

Now the key thing to notice is that a mode formed as a linear combination of phases 
g-iw Ar ^j^]-^ oil ^j' positive will be analytic and bounded in the lower half complex Ar 
plane. Conversely, if we take all the w' negative the linear combination will be analytic and 
bounded in the upper half complex Ar plane. Of course, neither such behavior is found for 
our mode s^j (which again, vanishes for positive Ar), indicating that we have both positive 
and negative frequency components in our mode. On the other hand, if we analytically 
extend the solution ( l22l) to negative x values by putting a branch cut in either the upper or 
lower half complex x plane, then we will have purely positive or negative free fall frequency 
solutions. At negative x, these solutions have the form 

i\T __ +tt4:Mui —iujT i4:MuiLog[—x] (2'^) 

(for positive x they of course agree with (l22l)). 

Our purely outgoing mode s^ can thus be decomposed as 

1 ^g.4M..p^_g-.4Af^^^]^ (24) 



g7r4MLj _ q-ttAMlo 



and the negative free fall frequency part is then 

-TTiMoJ 

^'^- = g-7r4Afa; _ g7r4A/<^ ^'^ " (2^) 

The complete wave-packet s is then decomposed mode by mode to obtain the negative free 
fall frequency piece s_. To find the norm of s_, we introduce one additional wave packet, 
-Srefiect: -Sreflect is defined to be a reflection of s around x = at very early times. Thus at early 
times it is bunched up against the inside of the horizon, and vanishes outside. This wave- 
packet corresponds to a particle which forms near the horizon, but then falls inwards into 
the singularity. We similarly reflect the individual modes fl22]) . It is easy to check from the 
form ([9]) of the inner product that the norm of Srefiect satisfies < ^reflect, ^reflect >= — < s,s >, 
with an analogous expression for the modes. It then follows from these definitions that 

N^ = S^ + e+''^^^S,eflect., (26) 

which gives < N^, N^ >=< s^^, s^; > (1 - e''^*^'^). We then use ([25]) to obtain 

<'~^'->= l_l8M. <'^'>- (27) 



Using (HM and normalizing our wave-packets, this corresponds to a thermal distribution of 
outgoing particles at the temperature T = l/SvrM. 

III. GHOST CONDENSATION 



y 



The ghost condensate [7|] may be thought of as a degenerate case of an irrotational perfect 
fluid. Indeed, if we consider an action of the form 



S, = J v^P(X), (28) 

with X = V^ipV^ip for some scalar field ip, then the physical system which results is precisely 
equivalent to that of an irrotational perfect fluid with pressure P and density 2P'X — P. 
Here P is any generic function of X of our choosing. 
The action (|28|) yields an equation of motion 



2P'{X)Dip + 4P"(X) V^f/^V^f/^V^V^yp = 0, (29) 

and an energy-momentum tensor 

j.f,u ^ 2P'{X)V^ifV''^ - g'^'PiX). (30) 

A ghost condensate background is defined to be a constant X solution having P' = 0. By 
an appropriate choice of cosmological constant, we see from (!30l) that such solutions may 
be taken to have vanishing energy-momentum. In flat space they correspond to Lorentz 
violating backgrounds which we may take to have the form ip = A^t, and if we expand the 
fleld in small fluctuations around such backgrounds, 

(p = AH + tt, (31) 

the equation of motion at leading order becomes 

7f = 0. (32) 

Note that the vanishing of P' implies that the leading term in the energy-momentum tensor 
will be linear in the vr fluctuations, and as a result the theory violates the null energy 
condition. Now, generically we should expect there to exist higher derivative operators in 
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the action, and these can stabihze the equation of motion (l32l) . We could for example have 
a higher derivative term of the form 

•^2 = - / V-9 j-2 , (33) 

leading in flat space to the wave equation 

A'P"{A')n + a^ = 0, (34) 

where A is the UV cutoff. 

At low energies, it follows from (|34l) that the spatial derivative terms dominate over the 
time derivative terms, and we have 

A'P"{A')n + a^ = 0. (35) 

In fact, this low energy form of the wave equation in fiat space is generic and does not 
depend on a specific assumption for the higher derivative operators such as ( !33l) . 
To keep expressions simpler we will henceforth adopt a specific form for P{X): 

P{X) = ^,{X-Ar. (36) 

This choice does not influence the results since we are interested only in the general form of 
the wave equation. The ghost condensate background in fiat space is then 

(p = Ah, (37) 

and the low energy dispersion relation following from ( l35l) is 

-' = «^- (38) 

We next consider the ghost condensate theory in the presence of a Schwarzschild black hole 
in 3+1 dimensions. The metric again takes the form (3), but now with the addition of the 
usual angular coordinates: 

ds^ = (1 - v^) dr^ + 2v drd r - dr^ -r^dQ\ (39) 

In the presence of the black hole, solutions of constant X having P' = continue to solve 
the equation of motion fl29|) . In this case V^ip now follows the world-lines of observers who 

11 



freely fall into the black hole from rest at infinity 13|. This corresponds to a background 
of the form 

^B = A'r. (40) 

Note that unlike in fiat space, the higher derivative terms in the equation of motion will 
perturb this solution and in fact lead to a gradual accretion of ghost condensate energy into 
the black hole. These effects are small however if we take 

1 A 

« -^ (41) 



2M ^ 

and 

A < Mpi, (42) 

which we will assume throughout the remainder of the paper. The first inequality is 
just the statement that the Schwarzschild radius of the black hole is much larger than the 
short distance cutoff of the effective field theory, which is certainly required for consistency. 
Having A much smaller than the Planck scale suppresses the Jeans instability in the ghost 
condensate, and also suppresses mixing between the vr excitations and the graviton. 

We may now expand the equation of motion (p9l) about the background (140|) to linear 
order in the it fluctuations. We obtain 

V^re^V^vr + erV/. V.vr + 2e^V^r V.vr = 0, (43) 

where we have defined C,^ = -^V^ipB to simplify expressions. As noted above, C,'^ follows 
the 4-velocities of observers falling freely into the black hole along geodesies, and as a result 
the third term in ( H3l) vanishes by the geodesic equation. This leaves us with 

V^reV^TT + eC^>.^urt = 0. (44) 

This wave equation follows from an action for vr of the form 

5r = [ v^g^e^.^f- (45) 



Now, we also have to include a higher derivative term in this action in order to stabilize 
the fluctuations. In the next section we turn to a discussion of the choice of this higher 
derivative term. 
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IV. HIGHER DERIVATIVE OPERATORS 



As pointed out in section (]III|) the choice of higher derivative operators typically does not 
affect the qualitative features of the low energy effective field theory of the ghost condensate 
particles. Unfortunately, just about all calculations performed to date of Hawking spectra, 
including the one reviewed in section ([TT]), make references to arbitrarily short distance 
physics. The detailed calculation presented here will also have this unpleasant feature, and 
we will be having to choose some particular higher derivative terms to include in our action. 

Fortunately, there is good reason to think that the results of the usual calculations, as 



well as this one, are independent of most of the assumptions made in the UV [ll|]. In 
particular, one can imagine a calculation scheme to make the UV independence manifest, in 
the following way: 

Consider foliating Schwarzschild space-time into a family of space-like hyper-surfaces- 
"nice slices"- which have the property that their extrinsic curvature is never much greater 
than ~ 1/M. We also assume that at infinity the time directions normal to the nice slices 
agree with the time-like killing vectors of the asymptotically flat space-time. That a family 
of nice slices actually exists was demonstrated in |11]. 

Now, we can imagine following the evolution of the state of the black hole from before 
it was formed, to some very late time afterwards. Hawking radiation would appear as low 
energy outgoing particles at large distances from the hole. The time evolution in question 
will be generated by a time dependent Hamiltonian operator, since the time direction normal 
to the nice slices is generally not the killing direction. On the other hand, all length scales 
governing the time evolution are at least of order the Schwarzschild radius. It then follows 
from the adiabatic theorem that modes of energy much larger than ~ 1/M will not be 
excited if they started in their ground states. 

Now, in the case of a theory with spontaneous Lorentz violation we have to be careful, 
since the modes we have argued are unexcited have high energy as viewed by observers 
moving normal to the nice slices. What matters for the validity of the low energy effective 
theory of the ghost condensate, however, is that there be no modes excited which have large 
energy with respect to the direction picked out by V^ips- The key point, though, is that 
slices of constant (pB have all the required properties to be nice slices so long as we do not 
look deep inside the horizon. 
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The extrinsic curvature on a surface of constant ifB is given by' 
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K^^^''-^. (46) 

Here Dyj^ was evaluated using the Christoffel symbols for the metric (!39|) which may be 
found in appendix A. This shows that as long as one is only interested in processes occurring 
outside of the horizon at r = 2M, then surfaces of constant ifB have acceptable extrinsic 
curvature to serve as nice slices. That these slices also have the correct asymptotic behavior 
at infinity follows from expressions ([2]) and (l40l) . 

It thus follows that low energy effective field theory can be expected to be valid for the 
calculation of the Hawking quanta emitted by a black hole in the ghost condensate.^ For 
this argument we do however have to assume that whatever the true UV completion of the 
ghost condensate may be (if one exists), it does not allow signal propagation at speeds much 
greater than the speed of light. Similarly, for the purposes of our calculation using the vr 
wave equation, we will require that the higher derivative terms we add to the action leave 
us with a well defined light cone for signal propagation. 

Perhaps surprisingly, it is actually rather hard to find higher derivative terms which have 
this property. In flat space, if we write a wave equation as 

F{tdt,td,)7r = 0, (47) 

then it turns out that the requirement for strictly subluminal propagation is for F{uj,p)~^ 



to be analytic in the region Imcj > | Imp| iJ]. A simple class of higher derivative terms 



which leave the vr wave equation with the appropriate analytic property may be added to 
(H5l) to yield the family'^ 

where a and c are free parameters which we will assume to be of order one. 



^ See Wald "General Relativity" page 230 for a definition. 

^ Though not perhaps for determining the quantum state after evaporation as was argued for in the standard 

case using nice slices. 
^ Note that we expect the Hawking spectrum to be well approximated by a calculation in the free field theory 

of the ghost phonons so long as interactions may be neglected. This is valid so long as energies of excited 

modes are much below the cutoff, implying that evolution due to the nonrenormalizable interactions is 

suppressed. 
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In flat space (HHI) yields the dispersion relation 



2 



UJ 



2 " /,2 ^ \2 



A2(^-^)- (49) 

Note that at energies far below A there are modes satisfying the usual ghost dispersion 
relation fl38|l . At energies far above A however, the higher derivative terms in the wave 
equation dominate, and signals propagate at the speed c.^ It is nice that the parameter c 
gives us one check on possible UV sensitivity of our results- in fact we will see that the value 
of c does not affect the Hawking spectrum in ghost phonons to within the approximations 
we employ. 

It should be pointed out that, while we expect the Hawking spectrum to be independent 
of the UV physics, there is an additional subtlety: Although in fiat space-time, the low 
energy wave equation satisfied by the vr's does not depend on the choice of higher derivative 
terms, this is not generally true in curved space-time. At this stage it is conceivable that 
different higher derivative terms, leading to different low energy physics in curved space-time 
could lead to different results for the Hawking spectrum. In the next section however, we 
will give some general arguments that this is not the case. 

Finally we note that we have dropped a source term in the wave equation which arises due 
to the fact that in the presence of the higher derivative operators the background ( HOl) is not 
an exact solution to the equation of motion. The true background is a small perturbation 
to this one under the approximation fHT]) that the black hole is large, and thus the wave 
equation we are using is also valid in this regime. 

V. THE HAWKING SPECTRUM OF GHOST PHONONS 

Equipped with the action pHj) . and the metric fl39|) . we now set out to calculate the 
Hawking spectrum in vr particles, following the same line of argument presented for the 2D 
massless scalar in section ([TTl)- For simplicity we will restrict ourselves to calculating the 
flux in s-wave modes and so will ignore tt dependence on the angular variables. 



There are however negative energy states in the spectrum which we will have to treat carefully. 
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In the (r, r) coordinate system the wave equation following from (HHj) takes the form 






V 



;i - ^)d'r 



,3v 2v 

- - —)9r 



c^r' 



a 



,3v 



,3v^ 5v^ 

V 



2c^r 



)dr 



A2 Vc2^^ C^ ^ '' 



7v^ 9v^ 
+ 



2r 



2c^r 






T?(3 - l)(-)(^^ + ^9r + ^) -.d'r + -,drdr - (1 - -j)dt + 



A^'c"^ 



,3f 2t> , ^ ,3f^ 



2r 



2f 



■)5.+ 



+ 



2c2^2 



TT 



0, 



(50) 



2c^r 2r^ 

where we have made the usual change of variables in 4D from tt to tt = -. Although this 
expression certainly looks like a mess, the reader will be glad to know that it is about to 
simplify significantly. 

As in section ([TT]) we are interested in solutions to this equation corresponding to modes 
which are purely outgoing at infinity and which vanish behind the horizon (note that de- 
pending on the value of c the sound horizon for the phonons actually is located at different 
radii, and it coincides with the usual black hole horizon for c = 1). Again we will look for 
solutions of the form 

^ = e-^'^^e'^^''^'''^'^'' (51) 



and hope that a WKB approximation will be valid (c.f. (TT8l) and (TT9|) ). The WKB approx- 
imation allows us to drop terms in fl50l) in which an r-derivative acts on v instead of on vf . 
Equivalently, since |^ ~ ^, we may also drop terms in which we have a factor of ^ instead 
of an r-derivative on tt. The fact that this approximation is indeed valid will follow once 
we have obtained an expression for k{r). Its validity will turn out to be equivalent to the 
requirement -^ <^ -^ that the black hole is much larger than the short distance cutoff. 
With the sub-dominant terms dropped we obtain the much simplified wave equation 



a 



d', + 2vdrdr + V'di + — -9^ + —dA 



2v 



2 \ 2 

^)9^ 



IT 



0. 



(52) 



A2 \c^ 
With the ansatz (15T1) and the second piece of the WKB approximation (iT9l) . this becomes 

{u;~vkf = ^^{e-^^{u;-vk)r- (53) 

Solving for k{r) we obtain: 



k{r) 



-{±1 + ^^UJ)V±' 



Va 2 



Va 



Veil 



±l + ^^^)2 + 4^(l-g)(±.. + ^f 
2^fl-^) 



(54) 
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where the prime on ±' indicates that it is separate from the other instances of ± (there are 
4 solutions to fl53|) ). Now note that u is the frequency of the mode as seen by an observer 
at infinity, so that we may take u; ^ A. In addition, having a mode with positive frequency 
uj and real wave number corresponds to a specific choice of the signs, so that the expression 
for k becomes 

k{r) = ^ ^ ^ ^. (55) 

From this expression it is a straightforward exercise to confirm that the WKB constraints 
(ITSl) and (fT9l) are indeed satisfied. 

In fact our expression for k{r) already reveals the physical reason for the absence of the 
Hawking effect in this theory. First of all, if we consider the limit of very small uj, then 
we can see that k becomes of order the cutoff. A, when v becomes of order 1 (again, this is 
assuming that c and a are also order 1). That is. A; ~ A for r of order but not necessarily 
close to the horizon radius. For larger u this will occur at still farther radii. This is to 
be contrasted with the standard case considered in section (JTll), for which such substantial 
blue shifting of the mode doesn't occur until it is located perturbatively close to the horizon. 
Along the same lines, let us choose an energy E satisfying A^ E ^ 1/M. Then it follows 
from flT6l) and fl55l) that cug r^ E when -^ ~ -^. Note now that this is a result which does 



2Af E ■ 

not require pushing the limits of the effective field theory since it only involves energies much 
smaller than A. We thus see that the outgoing modes have energies much larger than 1/M 
very far away from the horizon, while simultaneously satisfying the WKB approximation. It 
follows that the negative frequency components of these packets are negligible; by the WKB 
approximation they appear just like packets constructed out of ordinary plane waves at the 
given radius, and with the appropriate free fall frequency. Since these outgoing packets are 
then composed of only positive frequency components with energies much larger than 1/M, 
we may conclude that (to within the approximations employed) ghost condensate quanta 
are not radiated by these black holes. ^ Again, this argument is independent of assumptions 
about the behavior of the theory at high energies above the cutoff, and in addition, it is 
also independent of the low energy effects of the choice of higher dimension operators: In 
the WKB approximation, the wave-packets propagate locally as if they were in flat space. 



Making this arg ument more precise would involve going through an analysis with wave-packets such as 



that found in 



15| 
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and in flat space, the low energy behavior of the modes is generic, as explained in section 

m- 

On the other hand, the reader may be wondering if one could actually follow through 
with the same sort of near-horizon, infinite-blue shift argument given in section ([11]) for the 
2D massless scalar. We will now demonstrate that this is in fact possible, and that it yields 
a confirmation of our result. 

As in section ([TTl), the action (HHj) leads to a conserved inner product on complex solutions 
to the wave equation. This inner product is given by^'' 



< Si,S2 > 









- s 



^2) 



(56) 



On a surface of constant r the normal vector is n^ = ^^, and the expression simplifies to 

< Si,S2 > = 



T=const 



5trv.32 - ^(^rv.st) (ds2 + (^ - i)cfv^VpS2 
+ Y^si (^^,e^,Ds2 + (^ - 1)((^ - i)v.rre^v.v,.2 

+ ^rre'^V.V^V^Ss + V,CDs2)] - {si ^ S2) 



(57) 



We may use this inner product as before to associate annihilation operators with solutions to 
the wave equation of positive norm.^^ We are again interested in evaluating < a^(s)a(s) > 
for a purely outgoing wave-packet s, and will do so by evaluating the inner product defining 



^'^ The appropriate inner product may be determined as follows: Our action is a bilinear in ir with a sum 
of terms of the form / DattDbtt for some differential operators Da and Db- Now consider a new action 
given by replacing each of these terms by J DasIDbS2- This new action has a U{1) symmetry under 
which si and S2 transform with the same charge. Let J^ be the associated conserved current. Then the 
inner product is given by J^ n^^J^, with the overall coefficient chosen to enforce the relation P^ . with an 
appropriate normalization for the creation/annihilation operators. 

^^ We may still use expression p^ to define a{s), and the canonical commutation relations will still yield (fT^ . 
Note however that one must be careful to use the correct canonical commutation relations appropriate to 
a higher time-derivative field theory. See e.g. 16[. 
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a{s) at an early time when the wave packet is bunched up against the horizon (for c 7^ 1, 
against the sound horizon) . We can then decompose the wave-packet into a basis of high 
free fall frequency plane waves, enabling us to apply the condition that our state contains 
no high energy excitations. 

It is important to notice, however, that due to the fact that our action contains higher 
derivatives with respect to time, there is a complication. At high energies above A there will 
be ghost modes present, and, being negative energy states, these modes have the opposite 
relationship between their norm in the inner product (1571) and the sign of their free fall 
frequency. ^^ We note that the presence of these ghost modes in our calculation is not a 
problem since they occur only above the high scale A, and moreover we are now working with 
a free field theory- the negative energy states do not lead to an instability. On the other 
hand, when we decompose our wave-packet into pieces of positive and negative norm, it is no 
longer sufficient to decompose into pieces of positive and negative free fall frequency. The 
positive norm part of the wave-packet will now include a negative frequency ghost part. In 
fact, we will now show that the negative norm part of the wave-packet is entirely absent near 
the horizon- the negative frequency piece we would usually associate with particle creation 
is actually composed entirely of ghost modes and thus has positive norm. Likewise, the 
positive frequency piece is actually free of the ghost modes. It is in this way that we will 
confirm our earlier analysis which showed an absence of phonon Hawking radiation in this 
model. 

Near the sound horizon at r = ^, the full expression (15^ for the wavenumber (with 
appropriate sign choices) becomes at leading order 

fcx ~^^1 + ^^-, 58 



where x = r —. The expression for our outgoing mode outside the horizon then takes 

the form 



--V^('+^)^°^f"l. (59) 



Positive and negative free fall frequency modes can be obtained by analytically continuing 
in the upper or lower half complex x plane, as before. At negative x these modes take the 



^^ Unfortunately in order to use standard terminofogy, we are now using the term "ghost" to refer to 



negative energy states, not to general perturbations of the ghost condensate background. Please see 17 1 
for a discussion of the quantization of this theory. 
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form 

F^ = e "^ e e "^ (60) 

I _ 2A/A . .■ 2AfA T „„r ^1 

JV^ = g+-7^g-*-V^^Log[-x]^ (61) 

where we have now dropped the sub-leading terms in ^. 

Again, a purely outgoing wave-packet which vanishes behind the horizon can be formed 
as a linear combination of modes of these types. We will now show, however, that when we 
decompose P^^ and A^^^ into definite free fall frequency plane waves near the horizon, P^ will 
be composed purely of regular plane waves, while A^ will be composed purely of "ghost" 
plane waves. 

In flat space the dispersion relation (15^ becomes 

-' = :^(*= - ^)^- (62) 

This has four solutions: 

^-±\/?^- ">'' 

The ghost modes correspond to choosing the — sign inside the square root. At low energies 
these modes have imaginary k, but for tu > A, they become regular plane waves, despite 
having negative norm in the inner product (!56|) . If we take w ^ A, we can approximate 
the modes as having 

fc~±-±(±'-^). (64) 

Now, the modes we will expand our wave-packet in will be high frequency plane waves in 
the vicinity of the horizon as seen by a free falling observer, having the form 

where R is a spatial coordinate orthogonal to r (and the angular directions) near the horizon. 
From the metric fl39|) we can write i? ~ x -|- cr, so that the plane waves of interest become 

We will use this form for the plane waves to explicitly calculate the components of s. Before 
we do so however, we will first give an argument which will show the same results with less 
work. For simplicity we will here consider the case c = 1. 
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Along outgoing null rays with AR ~ At, the positive frequency plane waves (l65l) have 
leading behavior 

TTrg ~ e 2^S 
i-^Ar 



Here the letters r and / label right going or left going modes, while the letters n and g label 
normal or ghost type modes. Now note that it follows from the metric fl39l) that outgoing 
null rays near the horizon actually follow paths of the form x ~ ^4(1 + j^), where At and 
the constant A here are both taken much smaller than the Schwarzschild radius. Along such 
paths the positive frequency mode P^^, from ( I6OII . has the leading behavior P^^ ~ q^ve '^ j^ 
follows that, when we decompose P^ into the four types of modes in (167|) . the left going modes 
cannot be present, and neither can the right going ghost modes. Only a superposition of the 
right going normal type modes will yield an overall factor of e* v^ ^}^ We thus conclude 
that Pa; consists entirely of normal type modes. 

The interesting point now is that an analogous argument reveals that N^ is composed 
entirely of ghost type modes. The reason is that the plane waves with negative free fall 
frequency may be obtained from their positive frequency counterparts fl67j) by complex 
conjugation. In particular, the negative frequency, right going ghost modes will have the 
behavior e'^v^ ^ along outgoing null rays. As advertised, we then conclude that the 
negative frequency part of the outgoing mode s is composed entirely of ghosts, and thus 
has positive norm, in contrast to the standard case. It follows that the Hawking effect is 
absent, at least to within the approximations we have been employing. 

We now turn to an explicit calculation of the components of s using the expression (!57|) for 
the inner product. We will use the fact that ^ -C A, as well as the WKB approximation. 
Note that the WKB approximation is valid for the 5^, as well as for the plane wave modes 



^^ One could imagine cancelling an incorrect r dependence along a null ray by using combinations of plane 
waves with different values of uj. This possibility is excluded, however, since we may choose any null ray 
we like; such a cancellation could not be made to work across all outgoing null rays simultaneously. 
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(166|) at high frequencies. The inner product can then be written^^ 

/" _ a _ 1 

< s^, S2 > = 47ri / dr[s;^(9^ + vdr)s2 - jY^idr + vdr)sj{ — {dr + vdrf - dl)s2 

Jr>2m ^^ C C 



' r>2m 



+ J^^*u^(9- + ^9r){^{dr + vdrY ~ d',)^ - (s~ * ^ h)], (68) 



where as usual the tilde indicates that we have multiplied by r, and we only integrate outside 
the horizon since s^^ vanishes inside. We will take S2 to be one of the high frequency plane 
wave modes ( l66l) . In evaluating the integral we will make use of the identities 



[-Sdr + Vd,)-" - 9^)52 = ±'^CU252. (69) 

& \/a 



In these expressions, uj2 is the frequency of S2 appearing in expression (166|) . and the ±' 
coincides with that of expression fl63l) . In the first identity, we have dropped a sub-dominant 
term proportional to uj. The largest term that we can get from the integral (!68l) will be 
proportional to cjI, and we will drop all terms with weaker UJ2 dependence. This is an 
approximation which can be made arbitrarily accurate by computing the integral at earlier 
and earlier times; as the outgoing wave-packet is compressed against the horizon it becomes 
composed of arbitrarily high uj2 oscillations. A cutoff on uj2 for plane waves with appreciable 
overlap with s^^ would appear if we actually formed a wave-packet to localize the integral 
slightly away from the horizon. At earlier and earlier times, however, this cutoff becomes 
arbitrarily high. This observation allows us in particular to immediately discard the single 
derivative terms in ( l68l) . Similarly, we will keep only leading terms in x. Putting all of this 
together, the inner product can be written 



< So;, 



S2 >= 47ri^ / T'^^ls*j2 + (1 ±' l)cuj2S2drsJ + ic^ S2dls^*. (70) 



We will now integrate by parts to put all of the r derivatives on S2. We can then use 

drS2 = ik2S2 to obtain 



< S^,S2 >-- 



47r^ / {±'ujI + A;2(l ±' l)c^2 + klc')sZ*~S2. (71) 

^^C Jr>2m 



^^ More precisely, to justify the near horizon approximation we are employing we should use the full wave- 
packets formed out of the various modes. This can be done straightforwardly and does not affect the 
analysis. 
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With expression (1M|) for the possible values of k2, this becomes at leading order 

< s^, S2 >= 47r^ / {±'ujI ± (1 ±' 1)ujI + ujI)S:,*S2. (72) 

Ac Jr>2m 

It then follows that the order C(j| contribution to the inner product is non-zero only when 
S2 is either a positive frequency, right going, normal type mode, or a negative frequency, 
right going, ghost type mode. As argued above, this again allows us to reach the conclusion 
that, to within the approximations employed, phonon Hawking radiation does not occur in 
this theory. 

VI. DISCUSSION 

A few points deserve further discussion. First of all, although we have found no Hawking 
radiation in the limit that -g^ -^ 0, the reader may wonder if it is possible that a normal 
spectrum will appear at a higher order in -p^. But in fact this is not possible. For example, 
if the spectrum actually had a Boltzmann distribution going like e~^^^ with T oc 1/M, then 
we could keep uiM fixed as we take -g^ to zero, and the Hawking effect would not disappear. 
Our calculation is perfectly valid in this limit however, since even with tu = we obtain 
fluctuations around the background which satisfy the WKB approximation. Note that we 
can similarly rule out a spectrum such as e~^l'^ . The existence of the non-trivial u; = 
solutions is interesting in and of itself- they are indeed time independent, and one might 
wonder if they can be thought of as hairs for the ghost condensate black holes. This is 
actually not the case, since they have divergent energy momentum at the horizon as seen 
by freely falling observers. 

We have thus reached the conclusion that the ghost condensate phonons do not have 
a thermal spectrum of Hawking radiation, and in particular, the radiation vanishes in the 
limit 2^ ^ A ^ Mpi. To this author's knowledge, this is the first known example of a 
completely non-thermal Hawking spectrum. 

It seems likely that the ghost condensate is incompatible with quantum gravity, and it is 
interesting that signs of this fact show up in the Hawking spectrum, without requiring any 
particular couplings to other fields. 

There are a number of outstanding issues which would be worthwhile to address in future 
investigations. First of all, it would be nice to show explicitly that the GSL is violated in the 
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ghost condensate through direct apphcation of its null energy condition violating properties. 
In particular, one could try to send lumps of net negative energy into a black hole, shrinking 
the horizon area while reducing the entropy of the exterior. It would also certainly be very 
interesting to find theories other than the ghost condensate in which the Hawking spectrum 
is abnormal. Along the same lines, one might wonder whether it is possible to violate the 
GSL with a fiuid which is a small modification of the ghost condensate, perhaps even one 
which does not quite violate the null energy condition. Work is ongoing on these various 
issues. 
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APPENDIX A: CHRISTOFFEL SYMBOLS 

Here we list the Christoffel symbols for Schwarzschild space-time in the {r, r, ^, 0} coor- 
dinate system defined in the text. We have defined v{r) = — a/2 M/r and v' = dv/dr. All 
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symbols not listed are zero unless required by symmetry. 



TT 


= 


v\' 


pr 

TT 
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—vv' 


pr 
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v' 
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vr 




= 


vr sin^ 9 


Kr 


= 


— (1 — v'^)vv' 


Kr 
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-v\' 


pr 
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vv 


■pr 
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-(l-t;2)r 
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— (1 — f^)r sin^6' 


^ Or 


= 


1/r 


-p(f> 


= 


1/r 


p0 


= 


l/tan6' 


pe 


~ 


— sin Q cos Q 
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